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I.  INTRODUCTION 


Chirped  pulse  amplification  has  provided  access  to  intense,  few  cycle,  near  infrared  laser 
pulses  for  the  last  30  years  [1].  It  has  been  proposed  that  backward  Raman  amplification 
could  provide  similar  access  to  the  multipetawatt  [2]  or  even  exawatt  regime  [3,  4].  While 
experiments  have  yet  to  reach  the  multipetawatt  regime  [5-9],  it  is  possible  that  backward 
Raman  amplification  could  be  a  practical  source  of  long-wavelength  infrared  with  terawatt 
peak  powers. 

The  development  of  terawatt  power  pulses  in  the  long-wavelength  infrared  is  being  driven 
by  strong- field  science,  including  advanced  proton  acceleration  [10],  high  harmonic  genera¬ 
tion  [11],  mid-infrared  supercontinuum  generation  [12,  13],  and  nonlinear  optics  [14]. 

There  are  several  paths  to  high  power  long-wavelength  infrared;  optical  parametric  ampli¬ 
fication  (OPA),  difference-frequency  generation  (DFG),  optical  rectification  (OR),  or  lasing 
using  C02.  Frequency  down-conversion  from  the  near  infrared  to  long- wavelength  infrared  is 
inherently  inefficient.  In  principle,  the  Manley- Rowe  relations  limit  the  conversion  efficiency 
to  the  ratio  of  the  photon  energies  ^lwir/^nir,  about  10%.  In  practice,  the  conversion 
efficiency  is  an  order  of  magnitude  less.  For  example,  ZnGeP2  (ZNP)  optical  parametric 
amplifiers  become  inefficient  beyond  5  pm  [15,  16]  and  10  pm  [17]  due  to  absorption  [17]  and 
group- velocity  mismatch  [18].  Difference- frequency  generation,  in  combination  with  an 
optical  parametric  amplifier,  can  generate  tunable  broadband  pulses  with  10’s  to  100’s  of 
microjoules  [19,  20],  but  it  is  limited  by  phase-matching  and  optical  nonlinearities.  Optical 
rectification  in  GaAs  of  femtosecond  Tksapphire  pulses  can  generate  7  pm  to  20  pm  [18],  but 
is  limited  because  optical  rectification  relies  the  pump’s  spectral  wings  [19]. 

Existing  high  power  C02  amplifiers  can  create  picosecond  pulses  with  joules  of  energy 
[13,  21]  at  9.4  pm  and  10.6  pm.  However,  these  amplifiers  require  a  combination  of  large 
high  pressure  systems,  expensive  oxygen  isotopes,  and  a  high  power  seed  pulse  to  sufficiently 
broaden  the  C02  gain  spectrum  for  the  generation  of  picosecond  pulses.  Furthermore,  high 
gas  pressure  severely  limits  the  system’s  repetition  rate.  If  sufficient  broadening  is  not 
achieved,  these  systems  create  pulse  trains  instead  of  individual  pulses.  Attempts  have  been 
made  to  overcome  the  bandwidth  limitations  of  C02;  one  used  the  negative  group  velocity 
dispersion  in  GaAs  with  a  heatwave  to  compress  C02  laser  pulse  trains  [13],  and  the  other 
by  the  combination  of  self-chirping  and  a  conventional  dispersive  compressor  [22], 


3 


This  suggests  that  a  compact  10  pm  source  capable  of  terawatt  powers  and  picosecond 
durations  would  be  a  key  scientific  and  technological  development.  Plasma-based,  back¬ 
ward  Raman  amplification  can  provide  the  bandwidth  that  is  difficult  to  achieve  in  a  C02 
amplifier.  The  mechanism  is  self-phase-matched  and  optical  nonlinearities  in  plasmas  are 
weaker  than  in  crystals.  Additionally,  the  pumping  C02  laser  can  be  narrow-band,  which 
removes  the  need  for  operation  at  high  pressure,  and  in  turn,  improves  the  prospects  for 
high  repetition  rate  operation. 

In  backward  Raman  amplification,  a  backward  propagating  seed  pulse,  with  frequency 
and  wavenumber  (uq  ~  cu0  —  u>p,ki  ~  —  k0),  stimulates  the  coherent  backscattering  of  a 
forward  propagating  pump  pulse  (ui0,k0)  off  of  a  plasma  (Langmuir)  wave  (uip,kp  fa  2 k0). 
This  enhances  the  plasma  wave  and  further  drives  the  energy  from  the  pump  to  seed  waves 
leading  to  exponential  growth  in  the  seed  energy  until  pump  depletion  or  another  saturation 
mechanism  occurs  [23]. 

This  process  is  illustrated  in  FIG.  1.  The  pump  pulse  (blue)  is  moving  in  the  forward 
direction  (left-to-right)  while  the  seed  pulse  (green)  is  moving  backwards  (right-to-left).  A 
uniform  plasma  covers  the  entire  region.  Both,  the  pump  and  seed  arc  linearly  polarized.  The 
seed  pulse  is  injected  once  the  pump  pulse  is  overlapping  the  plasma,  as  shown  in  FIG.  1(a). 
The  initial  seed  pulse  has  a  full  width  at  half  maximum  (FWHM)  duration  of  3ps  and 
intensity  of  1  x  108  W cm-2.  The  seed  pulse  grows  exponentially  at  the  small-signal  gain  rate 
for  Raman  backscattering  y0  =  aoy'aWp/4  =  2.20(Ao[pm]/0[Wcm^2])1/,2(ne[cuf  _3])1/4Hz, 
where  ao  is  the  pump  normalized  vector  potential,  cuo  is  the  pump  frequency,  cop  is  the 
plasma  frequency,  Ao  is  the  pump  wave  length,  Iq  is  the  pump  intensity,  and  ne  is  the 
plasma  density  [23].  The  results  of  exponential  growth  can  be  see  in  FIG.  1(a)  as  the  seed 
intensity  increased  by  a  factor  of  30  since  entering  the  plasma.  This  corresponds  to  an 
approximate  gain  rate  of  7  ~  1/3. 7  ps.  During  the  exponential  growth  regime  the  seed 
pulse  duration  will  lengthen,  reaching  a  maximum  duration  of  5.2  ps  in  this  example.  This 
is  the  result  of  enhanced  backscatter  due  to  a  build  up  in  the  plasma  wave  at  the  tail  of 
the  seed  pulse.  Figure  1(b)  shows  when  the  seed  becomes  sufficiently  intense  to  deplete  the 
pump.  At  this  point,  the  seed  is  operating  in  the  superradiant  regime  [24],  The  leading  edge 
of  the  seed  is  backscattering  enough  of  the  pump  that  it  shadows  the  trailing  edge  of  the 
seed.  This  results  in  temporal  gain  compression  and  can  be  see  in  FIG.  1(b)  and  1(c),  with 
seed  durations  of  3.3  ps  and  1.5  ps,  respectively.  Figure  1(c)  shows  that  pump  depletion  has 
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continued.  The  seed  intensity  is  growing  roughly  linearly  with  propagation  distance  as  it 
sweeps  up  the  pump  energy. 


FIG.  1:  The  pump  (blue)  and  seed  (green)  intensity  at  three  times,  (a)  At  t  —  97ps,  the 
seed  has  entered  the  plasma  and  is  undergoing  exponential  growth,  (b)  At  t  =  110  ps,  the 
seed  has  become  sufficiently  intense  to  almost  fully  deplete  the  pump,  (c)  At  t  =  123  ps, 
the  seed  has  continued  depleting  the  pump  while  temporally  compressing.  The  seed  will 
continue  to  grow  after  the  final  frame  and  will  reach  a  final  seed  intensity  and  FWHM 
duration  of  54TWcm-2  and  1.3  ps.  The  plasma  density  is  1  x  1017cm-3. 

Researchers  at  Princeton  proposed  a  plasma-based  approach  for  pulse  compression  that 
takes  advantage  of  a  nonlinear  regime  of  amplification,  superradiance,  where  the  pulse  would 
shorten  as  it  grew  and  it’s  intensity  is  proportional  to  the  square  of  the  number  of  emitters 
(plasma  electrons)  [24],  The  formalism  for  this  process  was  originally  developed  for  growth 
from  noise  [25].  To  access  this  regime,  a  superradiant  pulse  must  be  intense  enough  u>b  >  u)p 
and  short  enough  r  ~  tt/cub,  where  ojq  is  the  pump  frequency,  cob  =  2u0y/a0ai  is  the  natural 
frequency  of  an  electron  in  a  plasma  wave  (also  called  the  bounce  frequency),  r  is  the  pulse 
duration,  a0,i  are  the  normalized  pump  and  seed  vector  potentials,  e  is  the  electron  charge, 
m  is  the  electron  mass,  and  c  is  the  speed  of  light  [23]. 

An  early  proposal  suggested  the  possibility  of  backward  Raman  amplification  for  gener¬ 
ating  exawatt  laser  pulses  [3]  which  would  open  access  to  new  fundamental  physical  regimes 
[26].  In  order  for  chirped  pulse  amplification  to  reach  higher  power,  excessively  large  grat¬ 
ings  are  required  to  keep  the  laser  fluence  below  the  grating  damage  threshold,  typically  on 
the  order  of  0.1  Jem-2.  Plasma-based  amplification  avoids  this  issue  by  not  using  gratings 
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to  stretch  the  pulse.  Additionally,  material  damage  thresholds  of  GW  cm-2  are  avoided  by 
the  use  of  a  plasma  instead  of  a  crystal  as  a  gain  medium. 

The  critical  power  for  relativistic  self- focusing  of  a  plasma  is  Pcr  ~  17((Uo/cap)2GW.  This 
is  the  power  at  which  plasma  nonlinearities,  and  instabilities,  begin  to  become  significant 
relative  to  linear  effects.  By  having  the  amplification  occur  faster  than  all  other  instabil¬ 
ities,  pulses  above  critical  power  could  be  generated.  The  development  of  the  pi-pulse 
model  dominates  the  theory  of  backward  Raman  amplification  [3] .  This  uses  a  three- wave 
model  to  describe  the  pump,  seed,  and  plasma  wave.  The  effects  of  modulation  instability 
(hlamentation),  wave  breaking,  and  forward  Raman  scattering  on  amplification  have  been 
discussed  at  length.  It  was  suggested  that  the  ideal  place  to  operate  is  near  wave  breaking 
and  with  the  highest  pump  energy  that  avoids  hlamentation.  The  pi-pulse  amplitude  and 
energy  should  grow  linearly  with  propagating  distance.  Simultaneously,  the  pulse  duration 
compresses  inversely  proportional  to  propagation  distance. 

Three  dimensional  particlc-in-cell  simulations  have  been  carried  out  which  show  regimes 
where  backward  Raman  amplification  is  feasible  [2].  For  example,  at  wavelengths  of  700  nm 
and  10  nm  there  are  parameters  that  yield  PW  powers. 

Experimentally,  the  observed  peak  output  powers  are  60  GW  [27].  Other  work  has  re¬ 
ported  efficiencies  of  6.4  %  [5]  for  a  pump  wavelength  of  800  nm,  seed  wavelength  of  890  nm, 
and  plasma  density  of  1.7  x  1019cm-3. 

Previous  work  [4]  mentioned  the  possibility  of  scaling  to  10  pm  wavelength  but  for  param¬ 
eters  where  the  initial  seed  intensity  was  equal  to  the  pump  intensity.  While  this  is  feasible 
for  800  nm,  it  is  not  for  10  pm.  This  makes  the  exponential  growth  regime  very  important 
to  the  final  seed  output. 


II.  DESIGN  CONSIDERATIONS 

High  power,  seed  sources  at  10  pm  are  lacking.  Therefore,  the  seed  must  make  use  of 
both  the  small  signal  gain  regime  and  depletion  regime  for  backward  Raman  amplification 
to  be  effective.  The  small  signal  gain  regime  is  required  in  order  to  get  the  seed  into  the 
depletion  regime  at  which  point  a  significant  fraction  of  the  pump  energy  can  be  transfered 
to  the  seed  pulse. 

The  pump  intensity  and  plasma  density  are  major  design  considerations  as  both  directly 
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influence  the  gain  rate,  wave  breaking,  and  other  instabilities.  Figure  2a  shows  a  contour 
plot  illustrating  various  design  limitations  in  this  plane.  The  gray  lines  show  the  time- 
independent  backward  Raman  gain  length,  the  length  over  which  the  seed  field  strength 
will  increase  by  a  factor  of  e1,  given  by  c/70,  where  y0  =  y/cJooy20/4  is  the  collisionless 
gain  rate.  The  black  lines  show  the  effect  of  electron-ion  collisions  on  the  gain  length  c/7, 
where  7  is  the  collisional  gain  rate  given  by  Eq.  (A4)  [23,  28].  The  electron- ion  collision 
rate  depends  on  the  plasma  temperature  and  ponderomotive  energy  with  a  functional  form 
approximated  by  uei  ~  (T  +  Up)~3^2 .  While  time- dependent  effects  are  important  [3],  the 
time- independent  gain  rate  is  illustrative  for  understanding  the  trade-offs  in  pump  intensity 
and  plasma  density. 


Xo=10p,m 


1012  1013  1014  1016  1016 

pump  intensity  (W/cm2) 


FIG.  2:  Contours  (black)  of  the  steady-state  Raman  backscatter  gain  length,  c/7,  are 
shown  for  2a  Ao  =  10  pm  and  2b  Ao  =  0.8  pm.  The  gain  length  in  a  collisionless  plasma  are 
shown  by  the  gray  lines.  The  red  region  shows  where  there  is  no  gain  due  to  collisions 
7q  <  vevs/A.  The  plasma  is  singly  ionized  with  temperature  T  =  leV,  In  A  =  10 ,  Z  =  1. 
The  blue  region  shows  where  wave  breaking  would  reduce  the  efficiency  once  the  seed  is  in 
the  depletion  regime.  The  green  region  shows  when  two-plasmon  decay  would  deplete  the 

pump. 

Shorter  gain  lengths  are  produced  with  higher  pump  intensities  and  plasma  densities 
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because  both  enhance  the  plasma  wave.  At  high  plasma  density,  electron-ion  collisions  {ye  oc 
ne )  will  damp  out  the  plasma  wave  faster  than  growth  from  backward  Raman  amplification 
(70  oc  Ue/4)  [23,  28].  The  difference  in  the  collisional  and  collisionless  gain  lengths  (black 
and  gray  curves  of  FIGS.  2a)  show  where  collisional  damping  becomes  important.  The 
red  lines  mark  the  instability  threshold  where  the  gain  is  zero  due  to  collisional  damping 
and  is  given  by  7q  =  ^ey,/4  (see  Appendix  A)  [23].  The  isocontours  of  the  collisional  gain 
length  (black  curve)  compress  against  the  boundary  where  collisions  dominate  gain  (red 
curve)  because  as  the  gain  approaches  zero,  the  gain  length  goes  to  infinity.  Additionally, 
the  collisional  damping  is  suppressed  relative  to  the  nominal  plasma  temperature  at  pump 
intensities  where  the  ponderomotive  energy  is  greater  than  the  plasma  temperature  J0  > 
1.07  x  1013  Wcm~2T[eV]/(A0[inn])2.  This  can  be  observed  in  FIG.  2a  when  the  red  line 
changes  slope. 

The  green  line  at  quarter  critical  density  ncr;t/4  fs  2.8  x  1020  cm_3/(A[pm])2  marks  where 
the  two-plasmon  decay  occurs  and  where  the  Raman  instability  becomes  absolute.  This  sets 
a  hard  upper  bound  on  the  possible  plasma  density  as  otherwise  the  pump  would  be  rapidly 
absorbed  by  the  plasma.  Typically,  other  instabilities  (forward  Raman,  parasitic  backward 
Raman,  or  hlamentation)  will  dominate  at  plasma  densities  below  quarter  critical  [2], 

Plasma  wave  breaking  can  limit  efficient  depletion  of  the  pump.  During  depletion,  each 
pump  photon  is  stimulated  to  scatter  into  an  additional  seed  photon  and  plasmon.  But 
the  plasma  wave  has  a  maximum  energy  density  that  it  can  support  [23],  in  other  words, 
a  maximum  plasmon  density.  Therefore,  when  the  maximum  plasmon  density  has  been 
reach,  the  seed  growth  is  stunted.  Plasma  wave  breaking  occurs  when,  locally,  the  seed 
has  encountered  a  pump  energy  equal  to  the  maximumly-supported  plasma  wave  energy 
of  that  region.  The  blue  curve  of  FIG.  2  marks  the  largest  pump  intensity  that  the 
plasma  can  support  during  100%  depletion  [29].  The  pump  intensity  should  be  below 
J0  =  4.7  x  10_15Ao[pm]  (ne[cm_3])3//2  Wcm-2  [29]  to  avoid  wave  breaking.  Plasma  wave 
breaking  does  not  affect  the  small  signal  growth  of  the  seed  pulse  but  the  overall  efficiency 
of  the  growth  when  the  seed  is  in  the  depletion  regime  [3].  Therefore  this  is  a  soft  limitation 
and  there  is  evidence  that  some  wave  breaking  is  preferable  as  it  can  limit  energy  growth  in 
the  tail  of  the  seed  pulse [3]. 

Within  these  constraints,  we  can  see  the  gain  length  varies  from  10mm  to  0.1mm.  There 
are,  however,  experimental  limitations  to  realizing  a  particular  plasma  length,  density,  and 
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uniformity.  In  particular,  available  sources  generate  long  wavelength  infrared  pulses  with 
microjoules  of  energy,  requiring  four  to  five  e-foldings  in  the  linear  growth  regime.  The 
need  of  using  both  the  linear  and  nonlinear  regimes  is  distinct  from  previous  studies,  which 
focused,  primarily,  on  near  infrared  or  visible  light  where  intense,  short  pulses  for  seeding 
are  already  available. 

The  parameter  space  for  10  pm  can  be  contrasted  with  FIG.  2b  which  shows  the  gain 
length  as  a  function  of  pump  intensity  and  plasma  density  for  800  nm.  It  is  clear  that  the 
operating  window  is  smaller  for  800  nm  but  there  are  not  the  same  technological  limitations 
on  seed  pulse  generation. 

In  summary,  this  motivates  using  laser  pump  intensities  of  around  1  x  1012  W  cm  2  and 
plasma  densities  around  1  x  101'  cm”3  in  order  to  avoid,  collisional  damping,  wave  breaking, 
and  have  a  gain  length  that  is  sufficiently  small. 


III.  NUMERICAL  RESULTS 

Motivated  by  the  regime  described  above,  one  and  two-dimensional  simulations  were 
carried  out  using  the  turboWAVE  framework  which  couples  a  finite-difference  time-domain 
electromagnetic  solver  with  a  collisional,  relativistic  cold  fluid  plasma  model  [30].  The 
simulation  domain  consists  of  several  sections  in  the  following  order;  a  “vacuum”  section, 
an  up  ramp,  a  uniform  plasma  section,  a  down  ramp,  and  a  final  “vacuum”  section.  The 
uniform  plasma  section  has  a  density  Hq  and  length  Lz.  The  “vacuum”  sections  have  an 
electron  density  10”4no-  The  “vacuum”  and  ramp  sections  are  all  1mm  long.  The  electron- 
ion  collision  model  requires  a  constant  plasma  temperature  and  is  set  at  1  eV. 

A  scaling  of  the  pump  intensity  for  fixed  plasma  density  was  carried  out  to  illustrate  the 
effects  of  collisional  damping,  the  exponential  gain  regime,  and  depletion  regime  on  the  final 
intensity  and  FWHM  duration  of  the  seed  pulse.  These  can  be  seen  in  FIG.  3.  For  this 
specific  set  of  simulations,  the  plasma  density  is  no  =  1  x  101'  cm”3,  the  plasma  frequency 
is  ojp  =  27r/(0.352ps),  the  grid  size  is  0.01/c”1,  and  the  time  step  is  0.009a;”1. 

The  pump  pulse  (Ao  =  10.6/;m)  enters  the  left  side  of  the  simulation  domain  at  t  —  0. 
The  initial  seed  pulse  (Ai  =  11.8/;m)  enters  the  right  side  of  the  constant  plasma  density 
region  as  the  pump  pulse  at  the  left  side  reaches  half  its  peak  intensity.  The  seed  pulse 
frequency  was  chosen  to  be  on  resonance  for  the  Raman  instability  oj\  =  ujq  —  u jp.  The 
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initial  seed  pulse’s  FWHM  duration  is  3  ps  with  300  pJ  cm-2  contained  within.  The  initial 
fluence  was  chosen  so  that  with  a  1  mm2  cross-section,  the  initial  seed  pulse  would  have  pJ 
energies.  The  plasma  length  is  Lz  =  10  mm.  The  pump  pulse  FWHM  duration  is  68  ps, 
which  corresponds  to  twice  the  plasma  length  2 Lz/c. 

Field  values  were  recorded  at  the  simulation  boundaries.  Spectral  box  filters  from  ±cup/2 
around  the  pump  and  seed  frequencies  were  used  to  extract  individual  field  envelopes  and 
intensity  profiles. 

Figure  3  shows  the  dependence  of  the  amplified  seed  on  pump  intensity  with  and  with¬ 
out  electron-ion  collision,  the  green  and  blue  curves  respectively.  At  intensities  below 
1  x  1011Wcm-2,  there  is  no  significant  gain  in  the  seed  pulse  because  the  gain  length  is 
comparable  to  the  plasma  length,  as  seen  in  FIG.  2a.  As  the  pump  intensity  is  increased 
from  1  x  1011Wcm-2  to  just  below  1  x  1012Wcm-2,  the  seed  intensity  in  the  collisionless 
plasma  grows  exponentially  with  pump  intensity.  This  means  that  an  increase  in  the  pump 
intensity  allows  the  seed  to  undergo  additional  e-foldings  within  the  fixed  plasma  length. 
This  stops  at  just  below  1  x  1012  W cm-2,  where  the  seed  begins  growing  more  slowly  with 
pump  intensity.  The  seed  pulse  has  grown  sufficiently  intense  that  it  is  beginning  to  deplete 
the  pump.  Similar  to  what  occurs  in  FIG.  1,  once  the  seed  pulse  begins  depleting  the  pump, 
the  rate  of  increase  in  the  seed  intensity  is  no  longer  exponential  but  roughly  linear  with  the 
encountered  pump  fluence.  For  fixed  plasma  length,  this  suggests  that  the  seed  will  grow 
roughly  linearly  with  pump  intensity.  The  rate  of  increase  in  FIG.  3  is  faster  than  linear, 
because  a  higher  pump  intensity  shortens  the  length  needed  for  the  seed  to  exponentiate  and 
reach  depletion.  This  increases  the  length  the  seed  can  deplete  and  hence  the  seed  grows 
faster  than  linear  with  pump  intensity.  For  a  pump  intensity  of  4  x  1012  W cm-2,  the  maxi¬ 
mum  seed  intensity  reached  5.5  x  1013  W  cnr2  with  a  duration  of  1.3  ps.  However,  we  were 
not  able  to  investigate  what  happens  as  pump  intensities  increase  above  4  x  1012  W  cm-2  be¬ 
cause  plasma  wave  breaking  causes  the  fluid  model  to  breakdown.  For  a  1  mm  cross-sectional 
area,  this  corresponds  to  a  peak  power  of  0.55  TW,  energy  of  0.72  J,  and  amplification  factor 
of  2.4  x  105. 

The  temporal  dynamics  of  the  seed  pulse  also  show  the  exponential  growth  and  depletion 
regimes.  When  the  gain  length  is  longer  than  the  plasma  length,  the  seed’s  final  duration  is 
equal  to  the  initial  duration,  as  seen  in  the  first  three  points  of  FIG.  3(b).  During  the  small 
signal  gain  regime,  the  point  of  maximal  growth  sweeps  backwards  at  c/2  [31],  essentially 
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FIG.  3:  The  final  seed  (a)  peak  intensity  and  (b)  FWHM  duration  are  shown  as  a  function 
of  pump  intensity.  The  blue  (green)  curves  show  the  result  of  neglecting  (including)  the 
effect  of  electron-ion  collisions.  The  maximum  seed  intensity  reached  is  5.5  x  1013  W  cm-2 
with  a  duration  of  1.3  ps  for  a  pump  intensity  of  4  x  1012Wcm-2. 


the  leading  edge  of  the  seed  is  driving  a  plasma  wave  so  that  increased  backscatter  will  be 
seen  later  in  the  seed  pulse.  This  can  be  seen  most  clearly  in  the  collisionlcss  plasma  where 
seed  durations  of  32  ps  are  observed.  When  collisions  are  included,  the  seed  duration  still 
grows  during  the  linearly  regime  but  only  at  higher  pump  intensities  where  the  electron-ion 
collision  rate  is  less  important.  Once  the  seed  begins  depleting  the  pump,  the  seed  duration 
rapidly  decreases.  What  is  occurring  is  not  compression  of  the  seed  energy,  but  amplification 
in  the  seed’s  leading  edge,  namely  temporal  gain  compression.  Pulse  durations  as  short  as 
1.3  ps  are  reached.  This  is  shorter  that  the  2ps-3ps  pulses  that  are  typically  created  in  high 
power  C02  lasers  [12,  22],  Previous  work  has  shown  that  some  degree  of  wave  breaking  is 
beneficial.  When  wave  breaking  occurs  after  the  peak  of  the  seed  pulse  has  passed,  it  can 
suppress  the  growth  in  trailing  pulses,  such  as  those  seen  in  FIG.  1(c)  [3]. 

A  scan  of  the  initial  pulse  duration  while  the  fluence  is  held  fixed  is  shown  in  FIG.  4. 
The  pump  intensity  is  2  x  1012  W  cm-2  and  all  other  parameters  are  the  same  as  FIG.  3. 
There  is  a  weak  dependence  of  the  seed’s  final  intensity  and  duration  on  it’s  initial  duration. 
This  is  promising  as  it  suggests  that  the  initial  seed  pulse  duration  is  not  a  critical  design 
parameter. 

The  intensity  maximum  in  FIG.  4(a)  can  be  understood  from  the  dynamics  of  the  small 
signal  gain.  The  point  of  maximal  gain  sweeps  backward  at  half  the  speed  of  light  [31].  If 
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Initial  Seed  Duration  (ps) 

FIG.  4:  The  final  seed  (a)  peak  intensity  and  (b)  FWHM  duration  are  shown  as  a  function 
of  the  initial  seed  duration.  These  are  fluid  turboWAVE  simulations  with  electron-ion 
collisions.  The  pump  intensity  is  2TW cm-2,  plasma  density  is  1  x  1017cm-3,  and  plasma 
length  is  10  mm.  The  black  lineouts  show  the  final  intensity  profiles  of  the  seed  pulse.  The 
lineout  axes  cover  a  temporal  and  intensity  range  of  80  ps  and  16  TW  cnr2. 

the  seed  it  too  short,  the  point  of  maximal  gain  will  pass  over  it  and  amplify  the  trailing 
tail.  If  the  seed  is  too  long,  the  head  of  the  seed  will  amplify  to  the  depletion  regime  before 
tail.  In  this  case,  the  pulse  would  have  reached  depletion  sooner  had  the  pulse  been  shorter. 
The  intensity  profiles  inset  in  FIG.  4(a)  show  this.  The  left  intensity  profile  shows  a  short 
pulse  on  top  of  the  amplified  and  temporally  lengthen  initial  seed  pulse. 

Two-dimensional  turboWAVE  simulation  results  are  shown  in  FIG.  5.  This  simulation 
is  similar  to  the  previous  parameter  scans  with  the  following  differences.  The  initial  seed 
intensity  and  duration  are  1  x  1010  W  cnr2  and  3  ps.  The  pump  and  seed’s  initial  e-1  field 
spot  sizes  are  0.4  mm.  The  pump  intensity  is  7.4  x  1011  W cm-2.  The  geometry  is  planar. 
The  longitudinal  and  transverse  spatial  coordinates  are  z  and  x,  respectively.  Along  the 
longitudinal  direction,  the  10  mm  plasma  is  constructed  the  same  as  the  one- dimensional 
simulations.  The  plasma  profile  in  the  transverse  direction  is  uniform  with  periodic  bound¬ 
ary  conditions.  The  simulation  domain  had  83328  and  64  cells  in  the  z-  and  x-directions, 
respectively.  The  cell  size  is  A z  =  0.01  k^1  and  Ax  =  2 k^1  in  the  z-  and  x-directions, 
respectively. 

The  depletion  of  the  pump  can  be  seen  in  FIG.  5(b).  The  seed  intensity  has  increased 
by  a  factor  of  360.  The  seed  has  a  smaller  spot  size  due  to  gain  focusing.  This  shows  that 
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FIG.  5:  Intensity  of  seed  (a)  and  pump  (b)  just  before  the  seed  exits  the  plasma  region. 


comparable  results  are  possible  in  the  two-dimensional  simulations  in  one.  The  transverse 
grid  size  is  too  large  to  resolve  plasma  perturbations  needed  to  drive  the  hlamentation 
instability.  Further  investigation  is  needed  to  determine  the  importance  of  that  effect. 


IV.  DISCUSSION 

Two  modeling  limitations  have  been  observed  when  simulating  backward  Raman  am¬ 
plification.  First,  the  fluid  simulations  are  limited  by  wave  breaking  of  the  plasma.  The 
most  promising  cases  for  amplification  tend  to  occur  when  the  pump  intensity  is  at  or  above 
the  wave  breaking  limit,  as  suggested  by  FIG.  2a  and  FIG.  2b.  Specifically,  this  occurred 
for  pump  intensities  of  4  x  1012  W cm-2.  The  pump  pulse  alone  will  not  drive  a  significant 
plasma  wave  by  itself.  As  the  seed  is  amplified,  it  beats  with  the  pump  and  drives  the 
plasma  wave  to  break.  Fundamentally,  this  limits  the  maximum  seed  intensities  that  are 
possible  to  reach  in  fluid  simulations  at  a  given  pump  intensity. 

Particle-in-ccll  (PIC)  simulations  offers  the  ability  to  model  the  amplification  process  in 
the  regime  above  plasma  wave  breaking.  Extensive  PIC  simulations  have  been  carried  out 
and  show  quantitative  agreement  with  the  collisionless  fluid  model.  But  when  the  number 
of  potential  e-foldings  is  sufficiently  large,  7 qLz/c  ~  4  to  5,  parasitic  Raman  backscatter 
competes  with  the  amplification  process  and  limits  growth.  The  seeding  of  the  parasitic 
process  is  several  orders  of  magnitude  larger  in  the  PIC  versus  fluid  simulations.  This  limits 
the  use  of  a  PIC  plasma  model  for  this  problem  because  it  is  of  the  practical  importance  to 
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starting  with  an  initially  weak  seed  pulse.  Detuning  the  interaction  with  a  spatially  varying 
plasma  density  [31]  may  suppress  amplification  of  noise  and  make  PIC  simulations  for  these 
parameters  more  feasible,  but  the  need  of  multiple  e-foldings  for  the  seed  to  reach  depletion 
makes  the  detuning  technique  more  challenging. 


V.  CONCLUSION 

The  turbo  WAVE  framework  has  been  used  to  carry  out  one  and  two  dimensional  FDTD 
electromagnetic  simulations  coupled  to  a  relativistic  cold  fluid  plasma  model  with  electron- 
ion  collisions.  Using  a  pump  pulse  that  could  be  generated  by  a  C02  laser,  it  was  shown  that 
a  seed  pulse  at  11.7  pm  could  be  amplified  to  5.4  x  1013Wcm-2  and  a  duration  of  1.3  ps. 
When  compared  to  the  initial  pump  pulse,  the  final  seed  pulse  is  10  times  more  intense 
and  50  times  shorter.  The  final  seed  amplification  is  weakly  dependent  on  the  initial  seed 
duration  which  is  promising  as  sources  are  limited.  The  ultimate  limits  on  efficiency  could 
not  be  determined  because  of  numerical  difficulties.  What  is  needed  is  a  collisional  kinetic 
model  in  which  the  noise  source  can  be  controlled. 


Appendix  A:  Temporal  Weak  Coupling  Gain 


As  derived  elsewhere  [23,  28],  the  Raman  backscattering  dispersion  relation  for  the  plasma 
density  perturbations  with  frequency  u j  and  wavenumber  k  is 


u2  —  u2k  +  iveu\ 


(cj  -  ujQf  -  [k  -  A'o  )  c2  -  ul  +  ivs  (u  -  uj0 ) 


uJ2k2v%s 


(Al) 


where  the  Bohm-Gross  frequency  is  cn2fc  =  oo2  +  3 k2v2h,  the  electron  thermal  velocity  is 
nth  =  \JT /m,  the  electron  temperature  is  T,  the  energy  damping  rate  of  the  electron  plasma 
wave  ue,  the  pump  frequency  and  wavenumber  are  cuo,k0,  the  energy  damping  rate  of  the 
electromagnetic  wave  (inverse-bremsstrahlung)  is  us  m  {ujp/ujQ)2ue,  and  the  quiver  velocity 
is  vos  =  eA0/(mc),  where  the  pump  vector  potential  is  AL  =  xA0cos(k0z  —  co0t)  [23]. 

The  scattered  electromagnetic  wave  should  have  a  frequency  u>s  ~  oj0  —  uip.  A  frequency 
detuning  of  Ac o  from  the  resonant  frequency  oj0  —  uek  defines  the  scattered  wave  frequency 
u js  =  uj0  —0Jek  +  Acu.  The  plasma  wave  should  be  approximately  the  Bohm-Gross  frequency 
except  we  will  define  it  with  a  real  frequency  shift  of  A u  for  detuning  and  a  complex  shift 
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of  Soj  for  gain  and  collisional  detuning,  that  is,  c o  =  u>ek  —  Ac u  +  8uj.  After  making  the 
approximation  that  8u,  Au  -C  coek  the  dispersion  relation  reduces  to 

(Su  +  ivs/ 2)  (5uj  +  ive/ 2  —  Au;)  =  — 7q,  (A2) 


where  the  resonant,  collisionless  gain  rate  is 

^2  =  Upk2v  os 

16cJe/c(cJo  ^ ek ) 

The  complex  frequency  shift  of  the  plasma  density  perturbation  is 


(A3) 


(A4) 


In  the  limit  of  no  collisions  ve,  vs  — y  0  and  resonant  Raman  Aui  — >  0,  Sco  =  i'y 0.  In  the 
limit  of  resonant  Raman  Aui  — y  0  but  no  Raman  gain  70  — >  0,  then  the  frequency  shift  just 
accounts  for  the  collisional  damping  of  the  plasma  wave  8uj  =  —ive/2.  In  the  limit  of  no 
Raman  gain  y0  — >  0  and  no  collisions  z/e,  vs  — >■  0,  then  there  is  no  induced  frequency  shift 
5u  =  0. 


If  the  interaction  is  resonant,  Au  =  0,  then  we  get  the  following  condition  for  an  insta¬ 
bility  7q  >  uevs/ A  [23].  If  ve  —  vs  —  0,  then  6co  =  Auj/2  +  i\J 7q  —  Acu2/4  which  says  that 
an  instability  only  exists  for  Acu  <  2yo.  The  power  spectrum  after  an  interaction  time  of  T, 
will  be  proportional  to  exp[2-\/7o  —  Acu2/4 T],  In  the  limit  that  the  detuning  is  much  less 
than  the  resonant,  collisionless  gain  rate,  the  FWHM  duration  of  the  seed  intensity  will  be 
approximately  tfwhm  =  ^/hi2T/70.  This  says  that  during  the  linear  growth  regime,  that 
the  pulse  durations  will  grow  with  time  due  to  gain  narrowing,  and  shorten  with  increased 
gain  rate  due  to  a  larger  gain  bandwidth. 


Appendix  B:  Collisional  Damping 

Collisional  damping  can  play  a  significant  role  on  backward  Raman  amplification  depend¬ 
ing  on  the  pump  intensity  and  plasma  density. 

The  intensity-dependent  electron-ion  collision  rate  is  given  by  vei  =  3z/0(uos/2uth)-3 
(3(uos/2uth),  where  the  zero  intensity  rate  is  u0  =  4v/27r-Z’2e4ne  In  A/(3m1/2T3'/2)  rs  2.91  x  1CT5 
?re[cm~3](T[eV])~3//2Hz,  the  quiver  velocity  is  vos  ~  25. 7  A  [pm]  y/7[ W  cm_2]cms_1,  the  ther¬ 
mal  velocity  is  uth  ~  4.19  x  10'  ^T[eV] cm s_1,  Q(x)  =  fQl  dzz2  [I0  (z2)  —  I\  ( z 2)]  e~z2 ,  and 
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the  functions  In(z)  are  modified  Bessel  functions  of  the  first  kind.  The  electron  charge  and 
mass  are  e  and  m.  The  Coulomb  logarithm  is  In  A  =  10  and  ionization  degree  is  Z  =  1.  The 
electron-ion  collision  rate  is  proportional  to  the  plasma  density  and  inversely  proportional 
to  the  electron  velocity  cubed.  The  thermal  and  quiver  velocity  play  a  role  in  the  overall 
collision  rate,  as  can  be  see  by  u0  and  3(uos/2uth)~3<5(uos/2r;th),  respectively.  This  has  been 
investigated  in  detail  [32],  As  laser  intensity  increases  and  the  quiver  velocity  becomes 
greater  than  the  thermal  velocity,  the  collision  rate  begins  to  decrease. 


l[W/cm2] 


FIG.  6:  The  electron-ion  collision  time  for  T  =  leV,  A  =  10  pm,  In  A  =  10  Z  —  1.  The 
green,  yellow,  and  red  regions  show  parameter  regimes  where  the  collision  time  is  good, 

mediocre,  and  poor. 


The  effect  of  laser  polarization  on  the  electron-ion  collision  rate  was  estimated  to  not 
cause  a  difference  larger  than  a  factor  of  3  in  the  rate  [32], 

Electron-neutral  collisions  are  not  considered  in  the  simulations  but  are  a  significant  con¬ 
sideration  if  the  plasma  is  not  fully  ionized.  The  collision  rate  is  given  by  ven  =  noaenvth 
where  uq  is  the  neutral  density,  aen  «5x  ICC15 cm2  [33].  To  include  the  increased  rate  of 
collisions  due  to  electron  quiver,  the  following  substitution  must  be  used  instead  of  the  ther¬ 
mal  velocity  alone  uth  — >■  y/u2h  +  u2s.  An  approximate  expression  for  the  electron  neutral 
collision  rate  is  uen  «  2.10  x  10_7no[cm_3]y/T[eV]  +3.76  x  10-13(A[pm])2/[Wcm-2].  Fig¬ 
ure  7,  shows  contours  of  constant  electron-neutral  collision  rate.  Generation  of  picosecond 
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duration  pulses,  requires  the  use  of  pump  intensities  and  neutral  densities  where  the  colli- 
sional  damping  is  slower  than  the  seed  pulse.  As  seen  in  FIG.  7,  this  means  neutral  densities 
below  1  x  101'  cm-3  would  be  feasible. 
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FIG.  7:  The  electron- neutral  collision  time  for  T  =  leV,  A  =  10  pm,  cross-section 
aen  =  5  x  10~15  cm2.  The  green,  yellow,  and  red  regions  show  parameter  regimes  where  the 
collision  time  is  good,  mediocre,  and  poor  relative  to  a  1  ps  pulse  duration. 


Finally,  previous  work  found  evidence  of  Landau  damping  [2]  but  it  saturated  quickly 
and  was  insignificant.  The  damping  rate  can  be  estimated  by  8uj/ujp  =  —  f  y^fcA^)-3 
exp  [—  A(/cA£>)-2]  [34],  where  the  characteristic  plasma  wave  wavenumber  is  k  —  2ui0/c 
and  the  Debye  length  is  Ad  =  \JT f  (47mee2).  For  a  pump  wavelength  of  A  =  10  pm, 
plasma  density  of  ne  =  lx  1017cm-3,  temperature  of  T  =  leV,  and  k\r>  =  9.33  x 
10 7 \JT [e V ]/ne[cm-3]/A [pm]  ~  0.03,  the  damping  rate  relative  to  the  plasma  frequency  is 
insignificant.  At  lower  densities  of  1  x  1015  cnr3  the  relative  growth  rate  is  Im [8uj\/ujp  ~  0.1 
and  could  warrant  further  consideration. 

The  units  in  this  section  are  cgs-Gaussian  unless  otherwise  stated. 
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